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Abstract 

We present a variational formulation of electrodynamics using de Rham even and odd differential 
forms. Our formulation relics on a variational principle more complete than the Hamilton principle 
and thus leads to field equations with external sources and permits the derivation of the constitutive 
relations. We interpret a domain in space-time as an odd de Rham 4-current. This permits a treatment 
of different types of boundary problems in an unified way. In particular we obtain a smooth transition 
to the infinitesimal version by using a current with a one point support. 
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f*^) ■ Introduction 

t^J- ■ A general framework for variational formulations of physical theories was presented in [1]. Appli- 

cations to statics and dynamics of mechanical systems appear in [2, 3]. In this paper we present a 
variational formulation of electrodynamics based on that framework. Our work is related to the gen- 
eral formulation of linear field theories in a symplectic framework contained in [4] and to the earlier 
formulations of electrodynamics contained in [5, 6]. Some of the results contained in this paper were 
k> ; announced in [7]. 

We are presenting a variational formulation of electrodynamics in an intrinsic, frame independent 
fashion in the affine Minkowski space-time using de Rham odd and even differential forms ([8, 9, 10]) 
which permit the rigorous formulations of electrodynamics and the description of the transformation 
properties of electromagnetic fields relative to reflections (see [5, 6]). Observable quantities like the 
charge contained in a compact volume or the flux of the electromagnetic field through a surface require 
the integration of differential forms. The use of odd quantities is not just a matter of elegance. Even 
if the space-time is assumed to be orientable, the alternative approach of using standard differential 
forms and the Hodge star operator require the choice of a specific orientation, i. e. the addition of an 
extrinsic structure. Other recent treatments of classical electrodynamics using odd and even differential 
forms can be found in [11, 12, 13]. 

Our construction is an example of application of the general variational framework [1]. Similar 
constructions are needed for the variational formulation of a general field theory. The linearity of 
classical electrodynamics and the choice of formulating it on the affine Minkowski space makes our 
presentation simpler. Relying on a variational principle more complete than the Hamilton principle 
our formulation leads to field equations with external sources. This variational principle also permits 
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the derivation of the constitutive relations which are usually postulated separately since the variations 
normally considered are not general enough to derive them from the variational principle. 

We interpret a domain in space-time as an odd de Rham 4-current. This permits a treatment of 
different types of boundary problems in an unified way. As an example we obtain a smooth transition to 
the infinitesimal version by using a current with a one point support. De Rham currents are essentially 
objects dual to differential forms. 

The present paper is organized in the following way. In the first part we provide the geometric 
structures needed for the rigorous formulation of electrodynamics. Part of this material is based on [5] 
and is briefly reported here for the sake of completeness. We recall the Cartan calculus for odd and 
even differential forms and their integration theory over odd and even de Rham currents. 

The second part contains the main results. We start with the construction of a suitable space of 
fields for electrodynamics (not a differential manifold) and a construction of tangent and cotangent 
vectors. A convenient representation of these objects is introduced. The definition of the space of 
fields is inspired by a similar construction suited for the statics of continuous media which is contained 
in the final section of [f ] . In Section 3 we formulate a variational principle for electrodynamics similar 
to the virtual action principle of analytical mechanics with external forces and boundary terms and 
derive the field equations which include the constitutive relations in addition to Maxwell's equations. 
The boundary problem in a finite domain is treated in Section 4. Section 5 contains the Lagrangian 
formulation of electrodynamics. The Lcgcndre transformation and the Hamiltonian formulation of 
electrodynamics in Section 6 conclude the paper. 

A. Preliminaries 

Here we provide the geometric structures needed for the rigorous formulation of electrodynamics in 
Part B. The material in Sections 1,2, and 4 is based on [5] and is briefly recalled here for the sake of 
completeness. Nevertheless, we add a more explicit presentation of some useful details. 

1. Orientations of vector spaces and vector subspaces. 

Let V be a vector space of dimension m ^ 0. We denote by F(V) the space of linear iso- 
morphisms from V to R m called frames. It is known that F(V) is a homogeneous space with re- 
spect to the natural group action of the general linear group GL(m,R) in F(V). Let GL + (m,R) 
and GL~(m,R) be the two connected components of the group GL(m,R). The set of orientations 
O(V) = F(V)/GL + (m,R) has two elements. This set is a homogeneous space for the quotient group 
H(m,R) = GL(m,R)/GL+(m,R). The sets E = GL+(m,R) and P = GL~(m,R) are the elements of 
the quotient group H(m, R) which is the group of permutations of the two elements of O(V). There is 
an ordered base (ei, e 2 , . . . , e m ) of V associated with each frame £ in a obvious way. 

Let W C V be a subspace of a vector space V. The subspace has the set 0(W) of orientations called 
inner orientations of W. Orientations of the quotient space V/W are called outer orientations of W. 
An outer orientation o" of W can be determined by specifying an inner orientation o of W together 
with an orientation d of V. Let (ei, . . . , e n ) be the base of W associated with a frame £ E o. This base 
can be completed to a base (e[, . . . , e' m ) of V with (e' l7 . . . , e' n ) = (ei, . . . , e„). The extended base can 
be chosen to be associated with a frame £' <E d . Let ir: V — > V/W be the canonical projection. The 
sequence 

(e'/,... ,e^_ n ) = (n(e' n+1 ),... ,tt(0) (I) 

is a base of V/W. It determines an orientation o" of V/W. Hence an outer orientation of W. The 
outer orientation o" of W constructed from o E 0(W) and d E 0(V) is the same as the orientation 
constructed from Po and Pd . 

We have introduced inner orientation of subspaces of dimension different from zero and outer orien- 
tation of subspaces of codimension different from zero. Integration theory of differential forms requires 
the possibility of assigning inner orientations to the subspace W — {0} C V and outer orientations to 
the subspace W = V. Two possible orientations are assigned to the subspace W = {0} C V one of 
which is distinguished. The distinguished orientation is denoted by (+) and the other orientation is 
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denoted by (— ). In agreement with the conventions established for orientation the outer orientations 
of the subspace W = {0} C V are the orientations of V. An outer orientation of W — {0} can be 
specified in terms of an inner orientation and an orientation of V. If the inner orientation is (+) and the 
orientation of V is o, then o is the outer orientation of W. The orientation Po is the outer orientation 
derived from (— ) and o. The subspace W — V has a distinguished outer orientation defined as an 
orientation of V/W. 

2. Multicovectors and multivectors. 

A q - covector in a vector space V is a mapping a: x g V x O(V) — > M. This mapping is g- linear and 
totally antisymmetric in its vector arguments. A q- covector a is said to be even, if 

a(v 1 ,v 2 , . . . ,v q ,Po) = a(v 1 ,v 2 , . . . ,v q ,o). (2) 

It is said to be odd, if 

a(vi,V2, ■ ■ ■ ,v q ,Po) = —a(vi,V2, ■ ■ ■ ,v q ,o). (3) 

The vector space of even g-covectors will be denoted by A q V* and the space of odd g-covectors will 
be denoted by A q V* . We will use the symbol A q V* to denote either of the spaces in constructions valid 
for both parities. The index p with the two possible values e and o will be used on other occasions. 

For the definition of exterior product of (even and odd) multicovectors we refer to [5]. Here we just 
recall that if two multicovectors a and a' are even or both are odd, the product a A a' is even. In other 
cases the product is odd. The exterior product is commutative in the graded sense and associative. 

Let {e K } K= i i ... ;TO be a base of V and let {e K } K= i ; ... >m be the dual base. Each element e K defines an 
even covector 

: V x O(V) — > K: (v, o) i — > (e K , v). (4) 

We choose an orientation o of V and introduce odd 0-covector e defined by e (o) = — e (Po) = 1 and 
the even covector e e defined by e e (o) = e e (Po) = 1. 

We come now to multivectors. We denote by K(x q V x O(V)) the vector space of formal linear 
combinations of sequences (vi,V2, ■ ■ ■ ,v g ,o) G x q V x O(V). In the space K(x q V x 0(Vj) we introduce 
subspaces 

A£00 = {Eti h{v\y 2 , o*) e K(xW x O(V)); 

E?=i^o(«i.«2."- .«a.o')=0 for each a G A q p V*} . (5) 

Subsequently we define quotient spaces A q V = K(x q V x 0{V)) / AP(V). Elements of spaces A q V and 
A q V are called even q - vectors and odd q - vectors respectively. A multivector is said to be simple if 
it is represented by a single element of the space V q x O(V) interpreted as a subspace of K(V q x 
0(y)). Evaluation of g-covectors on sequences (vi,v 2 ,--- ,v q ,o) G x q V x O(V) extends to linear 
combinations and their equivalence classes. If w is a g- vector represented by the linear combination 
J27=i ^i( v h v h ■ ■ ■ T v qT° l ) an( i a is a g- covector of the same parity as w, then 

(a,w) =E™=i A i«( v i. u 2>-- - .",,0*) (6) 

is the evaluation of a on w. We have constructed pairings ( , ): A q V* x A q V —> M. The exterior product 
of multivectors can be easily defined using representatives. The parity of the exterior product is odd 
if the parity of one of the factors is odd. It is even otherwise. The exterior product is commutative in 
the graded sense and associative. 

The left interior multiplications are the operations J : A q V x A q p , V* — > A q pp ^ q V* , defined for q ^ q' 
by (w J a, w') — (a, w A w'). The parity pp' which appears in this definition is constructed by assigning 
the numerical values +1 and — I to e and p respectively. The parity of the multivector w 1 must 
match the parity of the multicovector w J a. The right interior multiplications are the operations 
L : A q VxA q p W* -*■ A q pp , q 'v, defined for q > q 1 by (a',wL a) = (a'Aa,w). The parity of the multicovector 
a' in this definition must match the parity of the multivector w L a. 
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3. The Weyl isomorphism and a useful formula. 

The space A™V* is one-dimensional. This makes it possible to define the tensor product Af y<g>A™ V* 
as the set of equivalence classes of pairs (w, e) £ A q V x A" l V* . Pairs (u>, e) and (w' , e') are equivalent 
if there is a number A such that w' — Xw and e = Ae' or w = Aw' and e' = Ae. The equivalence class 
of a pair (w, e) will be denoted by w e. A tensor a e A'F <£> A™F* will always be presented as a 
product w®e. The set A q V® A" l V* is a vector space (see [5]). 

Proposition 1 ([5]). The linear mapping 

We q :A q e V ®A™V* -> A™- q V*:w®e^wAe (7) 

is an isomorphism. 

The mapping We g is called the Weyl isomorphism. 

We will show now that the values of any bilinear mapping b: A q V* x A% V* —> A™V* can be 
expressed using the exterior product. The following three technical lemmas will be used. 

Lemma 1. Let w e A q e V . Then 

w J (e D A e^; A ... A e" 1 ) = ^ (-l)Eli^-«)( e ^A...Ae^u I >e A<'+ 1 A...Ae^ (8) 

v 1 <...<v q 

where v q +\ < . . . < v m and {v q +i, • ■ • , v m ) denotes the complementary (m — q)-tuple of {y\, . . . , v q ). 
Proof: It is enough to prove the claim for a simple g-vector. If w = Wi A . . . A w q , then 

(u>i A ... A w q ) J (e Q A 6g A ... A e™) 

m 

= E(- 1 )' y " 1 « 1 ' u; i) w 9 J (••• J ( W 2J (e AeiA...Ae? A...Ae™))) 
i/ 1= i 

= ...= Yl (-l^-^sgnK... .^(e^) •••««,«,,) 

e G A e 1 A . . . A A . . . A e^ A . . . A e™ 

= E E (-l) E '= lK " <) sgna(e: < ' (1) ) w 1 )---(e^ ) )U , 9 ) 

e D A e 1 A . . . A A . . . A ee 9 A . . . A e™ 
= E E (-l^^^^sgna^'^O.-.^'^^eoAe^ 1 A...Ae^ 

;/!<.. .<v q a-eS(q) 

= J2 (-l)^- l(l/l ^ ) dct((e^, Ws ))e Ae^ 1 A... A<". (9) 
In the second line of this sequence of equalities we used the well known identity 

Q 

v J (a 1 A ... A a q ) = ^(-l)"" 1 ^, v)^ A...A?A...Aa« (10) 
i/=i 

where v is a vector and a 1 , . . . ,a q are covectors. In the third line we applied repeatedly this identity 
and we noted that when v\ < . . . < v q each missing covector e£* in the expression e D A e 1 A . . . A e v e x A . . . A 

ee q A ... A e™ occupied the {yi — i + l)-th place in the exterior product, otherwise if Vi-\-\ < v , < Vi-i, 
then it occupied the (yi — i — I + l)-th place, hence we get the factor (— l)S i= i( !/i ~^ sgn(^i, ... ,v q ) 
where the symbol sgn(^i, ... , v q ) denotes the sign of the permutation . . . , v q ). ■ 
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Lemma 2. If a e AjV* and w ® e G A«F ® A™V*, then 

(a, w)e = aA We g (w ® e). (11) 
Proof: It is enough to prove the claim for e = e Q A e\ A . . . A e™. By the Lemma 1 we get 
a A We ? (w ® e) 

= a A (to J (e Ae' A . . . A e™)) 

= ^ (-l)Eli("-- ! )( Wi ^i A . . . Ae^)a A e„ A e^ 1 A ... A e^™ 

ut<...<v q 

= E (-l)^^"'" ^-" 4 2 VV e Ae e 1 A...A^'Ae A<'+ 1 A...A<™. (12) 
vi<...<v q m<- ..<(u, 

The right-hand side of (12) reduces to 

£ (-l^L.K-^^-^a^.^e^e^ A...A<», (13) 
i/i < . . . < i/ ? 

since < . . . < v m and (f 9 +i, . . . , v m ) is the complementary (m — g)-tuple of (fi, . . . , Finally 
we note that moving each (for i = 1, . . . , q) to the i/j-th place in (12) requires — i transpositions, 
since v\ < ... < v q . Then each of the remaining Vi, i.e. those with i = q + 1, . . . , m, will necessarily 
be at the i^-th place. Hence we obtain 

dAWe,(s)«e)= S...,/"""'eoAee A -- A C ( 14 ) 

i/i < . . . < v q 



We denote by Hom(A™V*| Af V*) the space of linear mappings from A q e V* to A" l V* and by 

i 9 :Hom(A™y*| A« V*) -» A«V ® A"V* (15) 
the isomorphism characterized by 

(i,(0,o'(8»«) = (i(a / ),u> (16) 
for each Z G Hom(A™U*| A« V*), a' G A^F* and u G A™V. The pairing 

(, >:(A||V<8> A™V*) x (A«F* <g> A™V) -► R: (tu ® e, a <g> u) ^ (a,w)(e,u) (17) 

is used. 

Lemma 3. If I G Hom(A™V*\ A q e V*), then 

i(«)=aAWc,( l ,(i)), (18) 

for each a G A%V* . 

Proof: Let i q (l) — tuj ® e with w ; G A'V. Using Lemma 2 we obtain 

a A We q (i(l)) — a A We q (wi ® e) = {a, wi)e. (19) 

We will show now that 1(a) — (a,wi)e. Indeed if we denote by u G A" l V the dual basis of e then we 
have 1(a) = (1(a), u)e and 

(1(a), u) = (i(l),a ® u) = (wi ® e, a ® u) = (a, wi)(e,u) = (a, wi). (20) 
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We can now finally show that a useful expression can be obtained for the values of any bilinear 
mapping 

b:A«V* x A«V* -► A™V*. (21) 
We associate with b the linear mappings 

b:A q e V* -f A«V® A™V*:a~ v(b(o,-)) (22) 

and 

b: A«V* -► A«F ® A™F*: a' ^ i,(b(-, a'))- (23) 

Proposition 2. If b as in (21) is a bilinear mapping and b, b are the associated linear mappings (22) 
and (23), then 

b(a,a')=a'AWe,-(b(a))=aAWe,(E(a')). (24) 
for each (a, a') e A%V* x AgV*. 
Proof: Applying Lemma 3 to I — b(a, ■) we get 

b(a,a') =b(a,-)(a') = a' A We q > (i q > (b(a, ■))) = a' A We 9 /(b(a)). (25) 

The other equality is obtained in the same way by applying Lemma 3 to I = b(-, a'). ■ 



4. Integration of differential forms. Chains and currents. 

Let M be an afhne space modelled on a vector space V. A differential q-form on M is a differcntiable 
function A: M x X 9 V x 0(F) — > M depending on a point, q vectors and an orientation. It is g-linear 
and totally antisymmetric in its vector arguments. A differential form A is said to be even or odd if 
for each point in x G M it defines a even or odd multicovector, respectively. We note that a zero-form 
on M is a differentiable function /: M x 0(V) — ► R. The vector space of even differential q-forms will 
be denoted by $|(M) and space of odd differential q- forms will be denoted by <&£(M). The symbol 
$|(M) will be used to denote either of the two spaces when the distinction is of no importance. 

For the definition of the exterior product of two differential forms and for that of exterior differential 
of a differential form, we refer the reader to [5]. Here we just recall that if both forms A and A' are 
even or both are odd, the product AAA' is even. In other cases the product is odd. The parity of 
the differential dA of a g-form A is the same as the parity of the original form A. A g-form A can 
be interpreted as a q-covector field A: M — > A^V* . The exterior product and the exterior differential 
are extended to this alternative interpretation of forms. The left and right interior multiplications of 
even and odd multivector fields with even and odd multicovector fields are defined point by point in 
an obvious manner. 

A cell of dimension q or a q- cell in M is a pair (%, o), where x is a differentiable mapping \: W — > M 
and o is an orientation of V. For q = 0, MP is the vector space {0} with a single element 0. Hence a 
zero-cell in M is a pair of a point x £ M and an orientation of V. The integral of a g-form A on a 
cell (x, o) is the Riemann integral 

/ A= ■■■ A (x(si, • • • , s g ), Dix(si, • • • , s g ), ... , D g x(si, ... , s q ), o) dsi • • • ds q . (26) 
J(x,o) Jo Jo 

The integral of a 0-form / on a zero-cell (x, o) is the value J^ x „) / = /( x i <?)■ For each q we introduce 
the space K(X g (M)) of formal linear combinations of g-cclls. The formal linear combinations turn 
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into real linear combinations if cells are identified with elements of K(X q (M)). Integration of forms is 
extended to linear combinations by linearity. Subspaces IS|p(M) C K(X q (M)) are defined as the sets 

NP(M) = jc e K(X g (M)); J A = for each A £ $|(M)| . (27) 

Elements of the quotient spaces S£(M) = K(X q (M))/NP(M) are called even chains or odd chains of 
dimension q. We extend the sequence of even and odd chains to negative dimension q by defining the 
spaces E£(M) = {0} for each q < 0. A chain is said to be simple if it has a single cell as one of its 
representatives. Integrals of q- forms on q- chains are well defined. The integral of a g-form A on the 
class C oiC e K(X q (M)) is the integral of A on C. 

The boundary operator d assigns to a chain C £ S^(M) its boundary dC £ E^_ 1 (M). The boundary 
of a simple chain represented by a q- cell (x, o) is the chain represented by the combination 

£U(-1) < -V ,1) ,°)-(X (<,0) ,°)), (28) 
where the (q—1) -cells {x^' X \o) and (x^'°\o) are defined by 

x (*.*) : R9-i M:(si,... ,s,)i-»x(*i,-" ,*,) (29) 

with * replaced by 1 and 0, respectively. The cells introduced in (29) represent the faces of the simple 
chain. The construction of the boundary is extended to generic chains by linearity. The boundary of 
a boundary is the zero chain. It is known that Stokes theorem holds for chains and forms of the same 
parity. 

An even or odd de Rham current of dimension q on a manifold M is a linear function c: $|(M) — ► R. 
We will use the symbol J A to denote the value c(A). The spaces of forms are given certain topologies 
and the continuity of the function is required. Chains will be treated as currents. They form a dense 
subspace in the space of currents. We will consider only very simple examples of currents other than 
chains. Topological considerations are of little importance for these examples. The boundary of a 
current is defined by assuming that Stokes theorem holds for currents. Thus if c is a current of 
dimension q, then the boundary of c is the mapping 

0c:$« -1 (M)-fR:i4i-f / A=[dA. (30) 

Jdc Jc 

In addition to chains the odd de Rham current most frequently used is the Dirac current wS(x) of 
dimension m derived from an odd m- vector w and a point x £ M. If A is an odd to - form, then 

f A=(A{x),w). (31) 

J w5(x) 



B. Electrodynamics 
1. The space of fields. 

In this section we will construct a suitable space of fields for electrodynamics which has the same 
role of the space of motions in mechanics (see [3]). The space of fields is not a differentiable manifolds. 
Nevertheless, we will introduce in this space enough structure to permit the construction of tangent 
and cotangent vectors which are the objects needed in in the formulation of the variational principle. 
A class of admissible functions defined on the space of fields has also to be specified. This class needs 
to be large enough to include the action, which for electrodynamics is generated (in a sense that will be 
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specified) from a quadratic Lagrangian density, due to the linearity of the theory. Thus, only functions 
generated from quadratic mappings will be used in our variational formulation of electrodynamics. A 
larger class of admissible functions would be needed to deal with more general, non linear field theories, 
and the price of increased technical difficulties should be paid. 

Let M be the affinc Minkowski space-time of special relativity with the 4-dimcnsional model space 
V and the non degenerate metric tensor g : V — > V* of signature (1,3). The space of odd 4-currents 
with compact supports in M will be denoted by CM. Differential forms will always be presented as 
covector fields. 

We consider the set X($^(M); CM) of pairs (A, c), where c is an odd current of dimension 4 in M 
with a compact support Sup(c) and A is a local even 1-form 

A:U->AlV* (32) 

defined in an open set U C M containing the support of c. The 1 - form A will represent the electro- 
magnetic potential. 
A mapping 

k: M x A\V* x f\ 2 e V* -f A 4 Q V* (33) 
is said to be quadratic if for each x e M there exists a symmetric bilinear mapping 

S 2 k x : (a\V* x A 2 e V*) x (A\V* x A 2 e V*) - A 4 a V* (34) 

such that the mappings n x — n(x, •, •) and 5 2 n x satisfy the equality 

K x (a,f) = ^S 2 K x ((a,f)(a,fj), (35) 

for each (a, /) S A\V* x A 2 V* . We are using the standard definition of quadratic mappings in terms 
of polarizations. The polarization of the mapping k x is the mapping 8 k x and the equation (35) is the 
standard relation between a quadratic mapping defined on a vector space and its polarization. We will 
use the set of quadratic mappings (33) to introduce an equivalence relation in the set X($g(M); CM). 
In the next section they will be also used to define the admissible functions defined on the space of 
fields. 

Pairs (A, c) and (A',c') are equivalent if 

y,«° (x,A',dA') = J Ko(x,A,dA) (36) 

for each quadratic mapping n: M x A\V* x h 2 ^* — ► A 4 V* . The symbol x is used to indicate the 
identity mapping of M and also a point of M. 

If fx is an arbitrary odd 4-form on M and k is set to be the mapping 

k:M x AlV* x A 2 e V* -» A 4 Q V*: {x, a, f) ^ n(x), (37) 

then the equivalence condition (36) reduces to 

and implies that d = c. 

Equivalence classes of elements of X($J(M);CM) will be called fields. Our fields are similar to 
those used by Freed in [14]. The space of fields will be denoted by Q($g(M); CM) or simply Q. The 
equivalence class of (A, c) will be denoted by q(A, c) or simply q. There is a natural projection 

e: Q(<&*(M); CM) — > CM: q(A, c) c (39) 
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from the space of fields to the space CM of currents in M. It is not difficult to check that each fibre 
£ _1 (c) of the projection e is a vector space which will be denoted by the symbol Q($g(M); c) or Q c . 
Indeed, the sum of two pairs (A, c) and (A',c) with A: U -> A^V"* and A': [/' -> A*V* is the pair 
(A + A',c) defined in U H U' D Sup(c). Descending to the quotient with respect to the equivalence 
relation defined by (36) easily gives the sum in Q c . Therefore, the space of fields is the disjoint union 
of vector spaces Q = \J ceCM Q c - 



2. Functions, vertical tangent vectors and covectors in the space of fields. 

With each quadratic mapping n: M x h\V* x A 2 e V* — > A*V* we associate the function 

k: Q($g(M); CM) — > R: q(A, c) i-> ^ k o (a:, A, <L4). (40) 

The above mentioned admissible functions defined on the space of fields are those constructed in this 
way. They will be called differentiable. The space of such functions will be denoted by K($^(M); CM). 

We note that from the definition of the space of fields Q = Q($g(M); CM) (as a quotient space with 
respect to the equivalence condition (36)) follows easily that the differentiable functions separate points 
of Q, i.e. iik(q') = k(q) for each k G K(*J(M); CM), then q' = q. Indeed if k(q(A', c')) = k{q(A,c)) for 
each fc G K(®l(M); CM), then equation (36) holds for each quadratic mapping k: M x f\\V* x A^V* — > 
A^y*. It follows that c' = c and q(A', c) = q(A, c). 

The tangent space to a vector space Q c (i.e. to a fibre of e) coincides with Q c . 

The vertical tangent bundle to the vector bundle Q($^(M); CM) is the space 

VQ = Q x Q = {(q, 5q)eQx Q; e(q) = e(Sq)} (41) 

(e,e) 

where q — q(A, c) and 5q = q(SA, c) and A: U — > h\V* , 5 A: U f\\V* . The tangent projection is the 
canonical projection 

T Q :VQ^Q:(q,5q)~q. (42) 

There is no obvious choice of the bundle dual to \IQ. 

We will use the fibre derivatives of functions fc G K($g(M); CM) as models of covectors. The 
derivative Dfc of a function fc: Q — > R is defined for each current c separately. It is evaluated on a pair 
of vectors q = q(A, c) G Q c and 5q = q(5A, c) G Q c . The result is the expression 

Dfc {q{A, c),q{SA, c)) = ^fc (q(A + sM, c)) (0) 
as 



/* 9 

/ — k o (x, A + sSA, F + sSF) 
Jc ds 



(43) 

s=0 



where F = dA and 5F = dSA. The symbol dk(q(A, c)) will be used to denote the covector characterized 
by the pairing 

(dk(q(A,c)),q(5A,c))=Bk(q(A,c),q(SA,c)) (44) 

with vectors 5q = q{8A, c) G Q c . 

So we need to calculate the expression 



— n o (x, A + sSA, F + sSF) 

OS 



(x) = ~Dk x (A(x), F(x), 6A(x) © SF(x)), (45) 

s=0 



for each x G U. We are using the following known definition. If Vi,V~2,W are vector spaces and 
m: V\ x V2 — > W is a smooth mapping, then its derivative is the mapping 

Dm: V 1 xV 2 x (V x ® V 2 ) -» W: (a, /, <5a © Sf) ^ ^- (m(a + sSa, f + sSf)) (0). (46) 

ds 
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We define for each x E U the bilinear mappings 

A,: h\V* x A l e V* -+ A 4 a V*: (a, a') ^ 5 2 K x ((a, 0), (a', 0)), (47) 

» x :AlV* x AlV* - A^*: (a,/) ~ <5 2 ^((a, 0), (0, /)) (48) 

and 

*V A 2 e V* x A 2 e V* - A^y*: (/, /') ~ <5 2 ^((0, /), (0, /')), (49) 

obtaining the equality 

«x(o, /) = 2 A ^(°' a ) + M°> /) + 2 Vx (f' ^' ( 50 ) 

for each (a,/) £ A^y* x A 2 y*. The equality (50) in terms of the bilinear mappings (47), (48), and 
(49) will be useful to calculate the expression (45). 
We have the following lemmas. 

Lemma 4. If k x : A\V* x A 2 V* — ► A^V* is a quadratic mapping and \ x ,p x ,v x are the mappings 
defined above, then 

BK x {a,f,5a®Sf) = ^BX x {a,a,Sa® Sa) +Dfi x (a, f,Sa® Sf) + ^Du x {f,f,Sf®Sf), (51) 

for each (a,f,Sa®Sf) G A\V* x A*V* x a\V* © A 2 V* . 
PROOF: We have the equality 

1 1 

K x = -X x oAopri+fi x + -u x oAopr2, (52) 

where 

prv A\V* x A 2 e V* -» A^y*: (a, /) a, (53) 
A: AgV* — > AgF* x A\V*:a i— ► (a, a), (54) 
pr 2 : AlV* x A*V* - A 2 V*: (a, /) h- /, (55) 
A':A 2 e V* ^ A 2 e V* x A 2 e V*:a^ (a,a). (56) 
The claim follows from the equality (52) by noting that 

B Kx {a, /, Sa © Sf) = ^BX X o (A, DA) o (pn, Dpn)( 0) /, da © <f/) + D^(a, /, 5a © tf/) 

+ \pv x o (A', DA') o (pr 2 , Bpr 2 ){a, /, «5a © <5/) 
= ^DAa, o (A, DA) (a, Sa) + Du. x (a, f, 8a © Sf) 

+ ±Di/ x o(A , ,DA')(/,<f/) 
= ^DA x (a, a, <5a © Sa) + D/i x (a, /, <Sa © 5/) + ^Di/ X (/, /, Sf © 5/). (57) 

■ 

Lemma 5. IfV\,V 2 ,W arc vector spaces and 

b:V 1 xV 2 ^W (58) 

is a bilinear mapping, then 

D6(a, /, <5a © <$/) = b(a, Sf) + b(Sa, /), (59) 
for each (a, /, Sa © <$/) G Vi x V 2 x (Vi © V 2 ). 
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PROOF: From the definition (46) of the derivative it follows that 

D6(a, /, da ©<*/) = ^ (b(a, f + sSf) + sb(5a, f + s6f)) (0) 

- ^ (b(a, /) + sb(a, 6f) + sb(6a, /) + s 2 b(6a, Sf)) (0) 

= b(a,6f) + b(6a,f). (60) 



In view of Lemmas 4 and 5 the integrand in right hand side of the equality (43) reduces to 

Dk x {A(x),F(x),6A(x) ®5F(x)) = ^DX x {A(x),A(x),5A(x)®SA(x)) 

+ Dii x (A(x),F(x), 5A{x) © 5F{x)) 
+ ^Bu x {F(x), F{x),SF(x) © 8F(x)) 

= l - (X x (A(x), 6A(x)) + X x (SA(x), A(x))) 
+ ii x {A{x),5F{x)) + n x (6A(x),F(x)) 
+ i (u x (F(x),6F(x)) + v x {5F{x), F{x))) , (61) 



or 

Dk x (A(x), F(x),SA(x) © 6F(x)) = X x {A(x),SA(x)) + fi x (A(x), 5F(x)) 

+ fi x (SA(x),F(x)) + v x {F{x),5F{x)), (62) 

since \ x and v x are symmetric. By using the Proposition 2 contained in Section 6, this expression 
reduces to 

Bk x (A(x),F(x),SA(x) ®6F(x)) =6A(x) A Wei(A x (A(a;)) +6F(x) A We 2 (p x (A(x)) 

+ SA{x) AWei^(F(i)) +SF(x) AWe 2 {V x (F(x)) 
=SA(x) A (Wei(A x (A(aO) +%(F(x)))) 

+ SF(x) A (We 2 (-p x (A(x)) +V x (F(x)))) . (63) 

where the linear mappings X x and X x are obtained from the bilinear mapping X x as prescribed in 
Section 6, i.e., 

X x : t\\V* -» A*V © A>*: a » h(X x (a, •)), (64) 

J x : a\V* -» A*V © A^*:a h-> »i(A x (-, a)), (65) 

for each x £ U. The mappings Ti^., TZj, are obtained from [i x , and I7 X , 77 x are obtained from u x in the 
same way. 

We define from A^, X x the mappings 

A: M x h\V* -> A^V © A^V*: (a;, a) ^ A x (a), (66) 

I:M x A*V* -» A^8Ay*:(x,o) ^I x (a). (67) 
The mappings 71, 71, 7/, F are defined respectively from 7i x , ~p x , V x , V x in the same way. 
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By using the equalities (47) and (63) which hold for each x G U we obtain the equality 
d 



n -no(x,A + sSA, F + sSF) 

OS 



s=0 



= - (Wei ° (Xo(x,A)+JLo(x,F))) A 5 A 

+ (Wc 2 o (7t o (a;, A) + V o (x, F))) A SF. (68) 



s=0 



We have also used the graded commutativity of the exterior product. 
Hence, 

Dk(q(A,c),q(SA,c)) = f (x,A + s6A,F + sSF) 

Jc ^ s 

= - J (Wei o (A o (x, A) + p o (a;, F))) A <L4 

+ / (We 2 o (71 o (x, A) +Vo(x,F))) ASF. (69) 
Jc 

Since F = dA and SF = dSA, this equality can be converted to the form 

Dfc {q(A, c),q{SA, c)) = - J (Wei o(Ao(i,4)+po (x, dA)) 

+ d(We 2 o (71 o (x, A) + V o (ar,dA)))) A SA 

+ y d((Wc 2 o(7Io(x,A)+77odA))AM). (70) 

We note that the first integral contains the exterior product of SA and an odd 3-form, while the second 
integral contains the exterior differential of the product of 8 A and an odd 2-form. 

The obtained expression suggests the representation of covectors as equivalence classes of elements 
of the sets X($^(M) x $„(M); CM) introduced below. 

An clement of X(^ 2 Q (M) x & (M); CM) is a triple (G, J, c) of a local odd differential 2-form G: U 
A^*, a local odd differential 3-form J: £7 — > A^y* and a current c with support contained in J7. The 
odd 2-form G will be interpreted as the electromagnetic induction, while the odd 3-form J will be 
interpreted as the current. A covector p is an equivalence class p(G,J,c) of (G, J, c) e X(<&^(M) x 
$o(M); CM). The equivalence relations in X($£(M) x $^(M); CM) are based on the expression 



1 , „ . 1 

47TC~ 



— J ASA- — d(G AM) ) . (71) 



Elements (G, J, c) and (G', J', c') are equivalent if c = c' and 
for each (SA, c) e Q c . 

This equivalence relation is related to the equivalence relation (36). They can be used in the 
construction of dual objects. Indeed a (vertical) tangent vector might also (more generally) be defined 
as an equivalence class of pairs (SA, c) with respect to an equivalence relation similar to (36). A 
convenient representation of such relation is a relation similar to (72), but with (G, J, c) fixed, defining 
the equivalence between two pairs (SA, c) and (SA 1 , c). Such a construction is not needed in our case 
since each Q c is a vector space, so that its tangent space is canonically identified with Q c itself. 

The vector space II C of covectors associated with the current c will be used as the dual space of the 
space Q c thought of as the space of vertical tangent vectors with the pairing 

<P. Sq) c = JJ^JSSA- ^d (G A 5A)j , (73) 
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where 5q = q(6A, c) e Q c and p = p(G, J, c) G II C . The space of all covectors is the union 

n= (J n c . (74) 

cecm 

There is a natural projection 

e': II ^ CM: p(G, J, c)^c. (75) 
from the space of fields to the space CM of currents in M. The phase space is the space 

Ph = Q x n = {(q,p) e Q x n ; e(q) =s'(a)}. (76) 

(s,e') 

The symbol Ph c will denote the set Q c x n c C Ph. 



3. A virtual action principle for electrodynamics. 

In this section a variational principle for electrodynamics similar to the virtual action principle of 
analytical mechanics (see [3]) will be formulated. The construction is an example of similar construc- 
tions needed for the variational formulation of a general field theory. The linearity of the theory and 
the choice of formulating it on the affine Minkowski space makes our presentation simpler though 
containing the core of the general framework (which can be found in [1]). 

The action is the differcntiable function 

W: Q($g(M); CM) — > R: q(A, c) h-> J Lo(A,dA) (77) 
derived from the quadratic Lagrangian density 

L: AlV* x AtV* - A^V*: (a, f) ~ - JL</, f? e g-\f))y/\g{- (78) 

We are denoting by \f\g~\ the odd 4-covector defined by the Minkowski metric (see [5]). Later we 
will use the same symbol to denote the constant 4-covector field constructed from it. 

The 1-form A is called the potential, the 2-form F = dA is called the electromagnetic field. 

We note that the Lagrangian is a quadratic mapping which depends only on its second argument / 
and thus A = /j = and v does not depend on x in the formula (50). 

A phase ph = (q(A, c), p(G, J, c)) satisfies the virtual action principle if the equality 

(dW(q),6q)-{p,8q) c = (79) 

holds for each virtual displacement Sq = q(SA, c) G Q c . For each current c the dynamics associated 
with the current c is the set D c C Ph c of phases which satisfy the virtual action principle. The 
dynamics is the subset 

D= [J D c (80) 
ceCM 

of the phase space Ph. 

A phase space trajectory is a triple of local differential forms 

(A, G, J): U — > A\V* x A 2 D V* x A^V* . (81) 

The dynamics of a system can also be represented as a set D of phase space trajectories (A, G, J) such 
that for each current c with support included in U the phase ph = (q(A, c), p(G, J, c)) is in D c . 

The equation (79) is too abstract to be used directly. A more concrete expression is given in the 
following proposition. 
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Proposition 3. A phase ph = (q(A, c), p(G, J, c)) satisfies the virtual action principle if and only if 
the equality 

ir c J c ( d odA ) J ^\9\) ASA - d {{i^- 1 o dA) J y/\g\j A 5 A 



1 r rA 1 

-^JA6A- — 

C l 4*7TC 



2 ^M--d(GAM)] (82) 



is satisfied for each virtual displacement q(<5A, c). 

Proof: By applying the formula (70) to the action W we obtain that its variation is 

(dW(q),Sq) = J (-d(We 2 o V o dA) A 5A + d ((We 2 o V o dA) A 5 A)) , (83) 
since A = fi = 0. Moreover 

from which it follows after some calculation that 

V{x,f) = -^- c Alg- 1 {f)®^\, (85) 



for each / G A 2 e V* . Hence, 



\Ve 2 ovo (x,dA) = - JL ((A^g- 1 o dA) J vlffl) (86) 



and the variation of the action reduces to 



(dW(q), Sq) = ^ J (d ((Aig- 1 o dA) J v^|) A <5A - d (((A^g- 1 o dA) J y/\f\j ASA)) . (87) 
Recalling that, on the other hand 



(P,6q) c = jf (±JAM- JU(GAM)), 



(88) 



the claim follows. ■ 
A phase space trajectory belongs to the dynamics D, if and only if it satisfies the virtual action 
principle for each current c with support included in its domain of definition. There is a characterization 
of the dynamics of phase space trajectories in terms of differential equations. This is shown in the 
following propositions. 

Theorem 4. A phase space trajectory (A, G, J) belongs to the dynamics D if and only if it satisfies 
the Euler- Lagrange equation 

d((Alg- 1 odA) Jv /M) =^J (89) 

and the constitutive relation 

G=(Alg- 1 odA) jy/\f\. (90) 

Proof: If a phase space trajectory (A, G, J) satisfies the Euler-Lagrangc equation and the constitutive 
relation, then by substituting the expressions (89) and (90) of J and G in terms of the electromagnetic 
field F in the virtual action principle (82) it follows that (A, G, J) belongs to the dynamics D. 

The inverse implication will be proved in the next section. g 

The constitutive relation (90) produced by our variational principle corresponds to the momentum- 
velocity relation of analytical mechanics. 
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Proposition 5. A phase space trajectory (A, G, J) satisfies the Euler-Lagrange equation and the 
constitutive relation if and only if it satisfies the Maxwell's equations 



dG = — J (91) 

and the constitutive relation 

G = (Alg- 1 o F) J v^l, (92) 

with F = dA. 

Proof: The constitutive relation (91) is satisfied if and only if the equation (90) is satisfied for 
F = dA. 

If a phase space trajectory (A, G, J) satisfies the Maxwell's equations and the constitutive relation, 
then by substituting the expression (92) of the electromagnetic induction in the equation (91) we see 
that the Euler-Lagrange equation is satisfied, since F — dA. 

Conversely if the Euler-Lagrange equation and the constitutive relation are satisfied, then, again by 
substitution, we obtain that the Maxwell's equations (91) are satisfied. B 

We remark that since the virtual action principle (77) is more complete than the Hamilton principle, 
our formulation leads to the Maxwell's equations (91) with external sources. The proposed variational 
principle also permits the derivation of the constitutive relations (92) which are usually postulated 
separately since the variations normally considered are not general enough to derive them from the 
variational principle. 



4. The Dynamics in a compact domain. 

Let the current c consist in integrating an odd 4-form on a compact domain K C M with smooth 
boundary dK. Field configurations, tangent vectors and covectors are equivalence classes of equivalence 
relations based on the expressions 

no{x,A,dA) (93) 



IK 
and 

. \j A5A — d (G A SA) \ = 4r / JA5A-—[ G A 5 A. (94) 

K V c 47rc J c JK ^cJdK 

It follows that a field q = q(A,K) is represented by the restriction 

A\K:K->AlV* (95) 

of the potential A to the the domain K. A tangent vector Sq = q(SA, K) is represented by the 
restriction 

(5A)\K: K — > /\\V* (96) 

of the variation 5 A to the domain K. A covector p = p(G, J, K) is represented by the pair of the 
restriction 

G\dK:dK -► A 2 Q V* (97) 
of the electromagnetic induction G to the boundary dK of the domain K and the restriction 

J\K:k->hlV* (98) 

o 

of the current J to the interior K of the domain K . 

A phase is a pair (q,p) of a field q = q(A, K) and a covector p = p(G, J, K). 
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The action is 

W:Q K ^R:q(A } K) ^ f Lo(A,dA). (99) 
J K 

The virtual action principle is the equality 

(dW(q),Sq)-{p,Sq) K = (100) 

and the dynamics in the domain K is the set Dk C Ph of phases satisfying the virtual action principle. 

Proposition 6. A phase ph = (q(A, K), p(G, J, K)), defined in a compact domain K, belongs to the 
dynamics Dk if and only if the Euler-Lagrange equation 



d^A^odA) Jv^j) \K=^J\K (101) 

and the constitutive relation 

G\8K = ((A^g- 1 o dA) J V\g\) \9K (102) 

are satished. 

Proof: If q = q(A, K), p = p(G, J, K) and 8q = q(SA, K), then 



(dW(q),5q) = ±- J k d ((A^- 1 o dA) J y^f) A 6 A 



i- |^d((Alg- 1 odA) Jv^)aM 

" i /ax ( (A ' 3 " ° dA) J ^) A SA - (103) 



On the other hand 



{p, 5q) K = \ [ J ASA--?- [ G A SA. (104) 
c 2 7x Aire JdK 

Thus the virtual action principle assumes the form 

J k d ((A^g- 1 o dA) J v/M) A <5A - J ((A^g- 1 o dA) J ^g) A <5A 

= — / J ASA -[ G A 8 A. (105) 
By using variations with {SA)\dK = we derive the Euler-Lagrange equation 

/ \ O Ajr O 

d((Ae5 _1 °dA) J vlfll) |# = — JIA - . (106) 

Assuming that this equation is satisfied and using arbitrary variations, the constitutive relation 

G\dK = ((A^- 1 o dA) J v^f) (107) 

follows. ■ 
The following Proposition completes the proof of the Theorem 4. 
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Proposition 7. If a phase space trajectory 

(A,G, J):U -> a\V* x A 2 Q V* x A^V*. (108) 
belongs to the dynamics D, then the Euler-Lagrange equation 

d((Al 9 - 1 odA) j/M) =^J (109) 

and the constitutive relation 

G=(A 2 e g- 1 odA) j/M (110) 

are satisfied in U. 

Proof: If (A, G, J) is a phase space trajectory, defined in the open set U C M, which belongs to the 
dynamics D, then the equation 

d((Alg- 1 odA) A^\°K=^J\°K (111) 

and the boundary relation 

G\dK = ({Alg- 1 o dA) J v^|) \dK (112) 

are satisfied for each compact domain K C U. It follows that equations (109) and (110) are satisfied 
in every x S U. B 



5. The Lagrangian formulation of electrodynamics. 

The Lagrangian formulation of dynamics is the infinitesimal limit of the formulation in a compact 
domain with the domain shrinking to a point. A formal method which greatly simplifies the passage 
to the infinitesimal limit is to replace the compact domain — which is used exclusively as domain 
of integration — with the current c = 5(x)w, where S(x) is the Dirac delta function in x G M and 
w E AgV is an odd 4-vcctor, with ic^O. The construction of infinitesimal fields, tangent vectors and 
covectors is based on the expressions 

/ k o (x, A, dA) = (k (x, A(x), dA(x)) , w) (113) 

J S(x)w 



and 



8(x)w 



\ J A SA - ^-d (G A 5A)j = (^J{x) A SA(x) - -^—/■J- 6. lif.r). „■ ) . ( ] I J ) 



Since w ^ and A^V is one-dimensional, it follows from the first expression that a field q = q(A, c) is 
represented by the pair 

(A(x),F(x)) e A\V* x A 2 e V* (115) 

of an even 1-covector A(x) and an even 2-covector F(x). 
The second expression reduces to 

--!- ( f dG(x) - — J(x) ] A SA(x) + G(x) A 5F(x), w) , (116) 



47TC \ V C 

since d5A(x) — 5F(x). 
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It follows that a tangent vector 5q = q(5A, c) is represented by the pair 

(5A(x), 6F(x)) G a\V* x A 2 e V* (117) 
and a covcctor p — p(G, J, c) is represented by the pair 



G(x),dG{x) - — J(x) ) G A^y* x A^F*. (118) 



c 



The pairing 



assumes the form 



fa Sq) c = | c (^JAM- ^d (G A 5A)j (119) 



(P, *?) L = ( (dG(a;) - ^ A <L4(z) + G(x) A SF(x),w^ . (120) 

We have constructed the space of infinitesimal fields Q$ — A\V* x A^F* and the space of infinitesimal 
covectors Ws = A 2 V* x A^V* . Hence, the infinitesimal phase space is 

Ph s = Q$xIls = A\V* x A 2 e V* x A 2 V* x A 3 D V*. (121) 

The infinitesimal action is 

W(q(A, 6(x)w)) = (L(A(x),F(x)), w) . (122) 
The infinitesimal dynamics is the set 

Ds = {(q,Sq) G Ph s : V 5?eQ , (dW(q),5q) = {p,Sq) L } , (123) 

It is easy to verify that the infinitesimal dynamics Ds admits also the following more explicit expression 

Ds = {(a, /, .9, h) G Ph s , V (4M/)eAiv * xA j v * DL(a, /, 5a, */) = (h A oa + .g A */)} . (124) 

The infinitesimal dynamics Ds is characterized by the following Proposition. 

Proposition 8. An infinitesimal phase ph — (q(A, S(x)w), p(G, J, S(x)w)), with w ^ 0, belongs to 
the inhnitcsimal dynamics Ds if and only if the equations 



G(x) = (A 2 e g- 1 (F(x))) J/M (125) 

and 



dG(x) = —J(x) (126) 

c 



are satisfied. 



Proof: If q = q(A, 5(x)w)), p = p(G, J, S(x)w), and Sq = q(5A, 5(x)w), with w ^ 0, then the variation 
of the action (87), which can also be expressed in the form 

(dW(q),Sq) = -JL jf ((A^- 1 oF) J v^j ) A SF, (127) 

reduces to 

(dW(q),5q) = ~ (([Alg-\F{x))) J y^f) ASF(x),w) . (128) 
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Thus the virtual action principle 

(dW(q),5q} = {p,8qY (129) 

has the explicit form 

dG(x) - ^J(a0) A SA(x) + (G(x) - ((A^- 1 (F(a;))) J v4ffl)) ASF(x),w^ = 0. (130) 

Therefore if ph. satisfies the equations (125) and (126), then (130) is satisfied and ph 6 Dg. 
Conversely, if ph e Dg, then it satisfies the virtual action principle, which implies 

(dG{x) _ ^J(ao) A <L4(x) + (G(x) - ((A^" 1 ^*))) J y^f)) A <5F(z) = 0, (131) 

since w ^ and A^V is one-dimensional. The equations (125) and (126) follow, since 8A{x) and SF(x) 
are independent. g 

In a preceding section we showed that the dynamics of phase space trajectories can also be charac- 
terized by the Maxwell's equations and the constitutive relation. This fact can now be proved directly. 

Indeed, if a phase space trajectory (A, G, J) belongs to the dynamics D and a; is a point in the 
domain of definition of the trajectory, then the virtual action principle is satisfied for every infinitesimal 
current S(x)w. Thus the Maxwell's equations and the constitutive relation are satisfied in x, thanks 
to the previous proposition. It follows that they are satisfied in the whole domain of definition of the 
trajectory. 

To prove the inverse implication, we observe that the virtual action principle (82) can also be 
expressed in the form 

- JL I KV 1 °F)^\A6F = l^J,6A- (G A 6 A)) (132) 

or 

j {[p - (Alg- 1 ° F) J v4ffl) A SF + (dG - ^ jj A SA^j = 0, (133) 

for each virtual displacement q(5 A, c). Thus if a phase space trajectory (A, G, J) satisfies the Maxwell's 
equations, then it satisfies the virtual action principle for every current c with support contained in 
the domain of definition of the trajectory and hence it belongs to D. 




6. The Hamiltonian formulation of electrodynamics. 

We associate with the Lagrangian density 

L: a\V* x AlV* - A 4 a V*: (a, /) ~ -± (/, A 2 e g-\f)) VWl (134) 

the energy density 

E: AlV* x A 2 Q V* x A 2 e V* -» A A a V* (135) 

defined by 

E(a,g,f) = -^- c gAf-L(a,f) 

= -T L 5A/+-^(/,A^- 1 (.f))VH 



^ 5 a / + ^- / a(a^(/) J /m) 

±-(2g-Atg-\f)A^\g~\)Af (136) 
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and treat this mapping as a family 

E(a,g,-):A 2 e V* ^A 4 Q V* (137) 
of mappings on the fibres of the projection 

pr P : A\V* x A 2 Q V* x A 2 V* A\V* x A 2 a V* (138) 
onto the field-momentum space P = A\V* x A 2 V* . The set 

Cr(E,pr P ) = {(a,g,X) £ A\V* x A 2 q V* x A*V*; V 4A6A , v * DE(a, g, A, 0, 0, 5X) = o} (139) 

is the critical set of the family. The equality 

DE(a, g, A, Sa, Sg, SX) = — (Sg A A + g A SX) - DL(a, A, 5a, SX) 

= -^- c (5gAX + gA5X) + ^-8XA (a^" 1 (A) J 

= (<fo A A + (.g - (A^-^A) J ^|)) A SX) (140) 

implies 

D£;(a, ff ,A,0,0,JA) = -^( 5 - (^(Ajj^)) A <5A. (141) 
Thus we obtain the expression 

Cr(E,pr P ) = {(a, g, A) e A^* x A 2 a V* x A*V*; ff = A 2 ^ 1 (X) J ^\g~\} . (142) 
The critical set is the graph of the Legendre mapping 

A: A^y* x A 2 e V* - A^*: (a, A) ~ A^A) J VWl (143) 
For each a £ A\V* , the mapping A(a, •) is invertible. Its inverse is the mapping 

A(a, A^* - A 2 V*:g » A 2 e g (v^l L 5) , (144) 

where ^/|<? _1 | denotes the odd 4- vector characterized by = 1. It follows that the critical 

set is the image of the section 

o-.A\V* x A 2 a V* - A\V* x A^* x A 2 e V*:(a,g) ~ (a, 5, A^g (/F^ L .g)) (145) 

of the projection prp. The Hamiltonian density is the mapping 

£T = .Eoct:a£V*xa2V*->a2V*, (146) 

defined by the formula 

H{a,g) = -^ c gA( y A 2 e g( y ^/\g~ z ^\Lg)), (147) 

for each (a, g) £ A\V* x A 2 V* . The passage from the Lagrangian density L to the Hamiltonian density 
H is the Legendre transformation of electrodynamics. 
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We show that the energy density can be used to generate the infinitesimal dynamics D$. We consider 
the set 



lV*xfs 2 V*xA 2 V* 



De = {^^9,r) G Ph s ; 3 AeA , y , y {SMX)eKl 

DE(a, g, A, 5a, 5g, 5X) = — !— (r A 5a — f A 5g) \ . (148) 
This set is obtained by projecting the set 

^ = {(*>f>9,r,\) e Ph s x AlV*; ^^^^^ 

DE(a, g, A, 5a, 5.g, 5A) = — !— (r A 5a — f A 5g) \ (149) 

47TC J 

onto the phase space Ph 8 = A\V* x A^V* x A^V* Xj^JV*. 
It follows from (140) that X — f and that the set De reduces to 



D E = {(a,f,g,r,X)ePh 5 xA 2 e V*; X = f, 

V ( Sa ,S 9 ,S f) eKv* xa»v*xa 3 v* D ^< = ^ M 5« + .g A 5/)} . (150) 

It projects onto the infinitesimal dynamics 
D s =^{a,f,g,r)ePhs; V (4a , 



: 5/)eAiy*xA2y* 

DL(o, /, 5a, 5/) = - (r A 5a + g A 5f) J . (151) 

Hence, De — D$. 

It is clear from the definition of the set De that this set is included in the set 

{ (a, /, .9, r, A) G P/i 5 x A* V* ; (a, g, A) G CV(£7, pr P ) } . (152) 
The use of the mapping a results in 

De = {(a,f,g,r,X) e Ph s x A 2 e V*; X = f, 

V, s , * , * Dff(a, 5 ,5a,5g) = — (r A 5a - /A5 3 )l . (153) 

The Hamiltonian description of the dynamics 

£>s = {(«,/, 2, r)GPh 5 ; V (a0i4s)6Ail 



DII(a,g,6a,5g) = ^ (r A 5a - / A 5g)j 



|(o,/ )5) r)eP/i tf ; A^/^Ls) = /, r = o| . 



is obtained by projecting onto the phase space Phs- 
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(154) 



Proposition 9. An infinitesimal phase ph = (q(A, S(x)w), p(G, J, S(x)w)), with w ^ 0, belongs to 
the infinitesimal dynamics D$ if and only if the equations 

F(x) = Algo(y\j=t\\_G(xj) (155) 

and 



dG(x) = — J(x) (156) 
c 



are satisfied. 



Proof: We recall that an infinitesimal phase ph — (q(A,5(x)w),p(G,J,S(x)w)), with w ^ 0, is 
represented by 

A(x), F{x),dG(x) - —J(x), G(x)) (157) 



and use the Hamiltonian description (154) of infinitesimal dynamics Ds- The claim easily follows, g 
The resulting equations for the phase space trajectories (A, G, J) 

F = ^go(^\^\LG) (158) 

and 

dG = — J (159) 

c 

are called Hamilton's equations. The equations (159) are again Maxwell's equations and the equation 
(158) is the inverse of constitutive relation (92). 

7. References 

[1] W. M. Tulczyjew, The origin of variational principles, in the volume Classical and quantum inte- 
grability (Warsaw, 2001), 41-75, Banach Center Publ., 59, Polish Acad. Sci., Warsaw, 2003. 

[2] G. Marmo, W. M. Tulczyjew, P. Urbanski, Dynamics of autonomous systems with external forces, 
Acta Physica Polonica B, 33 (2002), 1181-1240. 

[3] A. De Nicola, W. M. Tulczyjew, A variational formulation of analytical mechanics in an affine 
space, Rep. Math. Phys., 58 (2006), 335-350. 

[4] W. M. Tulczyjew, A symplectic framework of linear field theories, Annali di Matematica pura ed 
applicata, 130 (1982), 177-195. 

[5] G. Marmo, E. Parasecoli, W. M. Tulczyjew, Space-time orientations and Maxwell's equations, Rep. 
Math. Phys. 56 (2005), 209-248. 

[6] G. Marmo, W. M. Tulczyjew, Time reflection and the dynamics of particles and antiparticles, Rep. 
Math. Phys. 58 (2006), 147-164. 

[7] A. De Nicola, W. M. Tulczyjew, A Note on a Variational Formulation of Electrodynamics, Proc. XV 
Int. Workshop on Gcom. and Phys., Tenerife (Spain), 2006, Publ. de la RSME 11 (2007), 316-323. 

[8] J. A. Schouten, Ricci calculus, Springer, Berlin, 1955. 

[9] J. A. Schouten, Tensor Analysis for Physicists, Oxford University Press, London, 1951. 

22 



[10] G. dc Rham, Varietes Differentiables, Hermann, Paris, 1955. 

[11] T. Frankel, The Geometry of Physics: An Introduction, Cambridge University Press, Cambridge, 
1997. 

[12] F. W. Hehl and Yu. N. Obukhov, Foundations of Classical Electrodynamics, Charge, Flux, and 
Metric, Birkhauser, Boston, MA, 2003. 

[13] I. V. Lindell, Differential Forms in Electromagnetics, IEEE Press, Piscataway, NJ, and Wiley- 
Interscience, 2004. 

[14] D. Freed, Classical field theory and Supersimmetry, IAS Park City Mathematics Series Vol. 11, 
IAS, Princeton, 2001. 



23 



